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Abstract For any positive integer n, let ¢(m) and S(n) be the Euler function and the Smaran- 


dache function respectively. In this paper, we use the properties and the curve figure of these 

two functions to study the solvability of the equation )*> S(i) = o( 2a), and prove that 
i=1 

this equation has only two positive integer solutions n = 1,10. 
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81. Introduction and result 


For any positive integer n, the famous F.Smarandache function S(n) is defined as the 
smallest positive integer m such that n divides m!. That is, S(n) = min{m: me N, n|m}, 
where N denotes the set of all positive integers. From the definition of S(n), it is easy to see 
that ifn = p{'p5?---p?* be the factorization of n into prime powers, then we have 


S(n) = mas ee Nhs 


It is clear that from this properties we can calculate the value of S(n), the first few values of 
S(n) are: S(1) = 1, $(2) = 2, $(3) = 3, S(4) = 4, S(5) = 5, S(6) = 3, S(7) = 7, S(8) = 4, 
S(9) = 6, S(10) = 5, ---. About the arithmetical properties of S(n), some authors had studied 
it, and obtained many interesting results. For example, Lu Yaming [2] studied the solvability 
of an equation involving the F.Smarandache function S(n), and proved that for any positive 
integer k > 2, the equation 

S(mz + m2 +--+ + my) = S(m1) + S(me2) +--+ + S(mz) 


has infinite group positive integer solutions (m1, me, --- , Mz). 
Jozsef Sandor [3] proved that for any positive integer k > 2, there exist infinite group 


positive integers (m1, m2, ---, mx) satisfying the inequality: 
Also, there exist infinite group positive integers (m1, m2, --- , mz) such that 


S(mi1 + M2 +--++ mg) < S(mz) + S(mz2) + +++ + S(me). 
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Rongji Chen [5] studied the solutions of an equation involving the F.Smarandache function 
S(n), and proved that for any fixed r € N with r > 3, the positive integer n is a solution of 


S(n)" + S(n)" 1 4---4+ S(n) =n 


if and only if 


where p is an odd prime satisfying p"~1 +p’? +---+1|(p—1)!. 

Xiaoyan Li and Yanrong Xue [6] proved that for any positive integer k, the equation 
S(n)? + S(n) = kn has infinite positive integer solutions, and each solution n has the form 
n = pny, where p = kn, — 1 is a prime. 

For any positive integer n, the Euler function ¢(n) is defined as the number of all positive 
integers not exceeding n, which are relatively prime to n. It is clear that ¢(n) is a multiplicative 
function. 

In this paper, we shall use the elementary method and compiler program to study the 
solvability of the equation: 


S(1) +5@) +--+ 5(n) = 6 (MED), (1) 


and give its all positive integer solutions. That is, we shall prove the following: 
Theorem. The equation 


S(1) + S(2)+---+S(n) =¢ (mae) 


2 


has and only has two positive integer solutions n = 1, 10. 


§2. Main lemmas 


In this section, we shall give two simple lemmas which are necessary in the proof of our 
Theorem. First we have the following: 
Lemma 1. For any positive integer n > 100, we have the inequality 


i 1? n? 
aay 
2, 5 $199 inn 
Proof. From the mean value formula of S(n) (See reference [7]) 
nw? a x 


we know that there exists one constant N > 0 such that 


n 2 2 2 2 2 2 
> S(i) 2 wT n | 1 7 n Z T n 
i=l 


S72 Inn 1199 12 Inn ~ 11.99 Inn 


holds for all positive integer n > N. We can take N = 100 by calculation. This completes the 


proof of Lemma 1. 
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Lemma 2. For Euler function $(n), we have the estimate 


1 +1 1 

3 (“e ') a mn ) el. = ae 

ln (2 In —a) 
Proof. Let n = p{'p$?---pp* be the factorization of n into prime powers, then there 
always exist some primes pi, po, ---ps such that pipo---ps; > n. From [1] we have 
x 
S"inp=2+0 (=). 

log x 


pK@ 


by this estimate we know that 


Inn< So inp; < S- Inp; < ps < 2Inn. 
i=1 PiX<Ps 
Thus 1 1 
So-< S- — <InInp, < InIn(2Inn). 
Pp Pi 


p|n Pi<Ds 


1 1 
Note that ¢(n) = n{[ (1 ), if mnt ) is even, then 
Pp 


pin 
n(n + 1) n(n +1) 1 
SNE ey _ nati) so 
o( 2 2 II 
p| 2a) 
In(1— 4)415_15 
= MHD) nh gga PTET 
4 
= » oer xu {In(1—4)+ +5] 
_— n(n bre 1) pein) 5 P j mF) 4 Pp Pp 
4 7 2 2 
> n(n +1) pi MED pga > 
4 
> Uo), ef. e7) Sin n(2In *S) 
4 
n(n+1) 3 1 
— -e4- 15 n(n+1) . 
: nen 


1 
If aint) is odd, we can also get the same result. This completes the proof of Lemma 1. 


§3. Proof of the theorem 


In this section, we shall complete the proof of our Theorem. First we study the tendency 
of the functional digraph 


KK pee) a 1 ng? 
a mn 5(2in 2@4D) 11.99 “Ina” 


By use of Mathematica compiler program we find that f(a) > 0, if a > 100754. 
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No. 


find 


programming language, we obtain that the equation (1) has no any other positive integer 


6 
7.5 x10 
sx1o® 


2.5 xo ® 


é 20000 40000 60000 60000 107000120000 
-2.5 x10 


-s>xa08 
6 
-2.5 x10 


figure 1 


From the figure 1 we know that if n > 100754, then 


” nm mn?  n(nt+i1) 3 1 n(n + 1) 
S(i) < oe . 
> OTe in a in (21m BED) . a( 2 ) 


If « € (100754, +00), we use Mathematica compiler program to compute f (x), then we 


that the derivative f’ (x) is positive, so (2) is also true if « > 100754. 


Now we consider the solution of (1) for all n € [1, 100754]. By use of the computer 


solutions except n = 1, n = 10. This completes the proof of Theorem. 
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The computing programme is given as follows if n € [1, 100754]. 
# include “stdio.h” 
# include “math.h” 
# define N 100754 
int S(int n) 
{int ret=1,num=n; 
unsigned long int nn=1; 
for(ret=1;ret<=n;ret++) { nn=nn«ret; 
if(an%num==0) break;} if (ret>n) ret=n; 
return ret; } 
int SumS(int n) 
{int ret=0,i; 
for(i=1l;i<=n;i++) ret+=S(i); 
return ret; } 
int coprime(int i,int n) 
{ int a=n,b=i; 
while(a!=b) { if(a==0) return b; 
if(b==0) return a; 
if(a>b) a=a%b; 


else 

b=b%a;} 

return a; } 

int Euler(int n) 

{int ret=1,i; 

for(i=2;i<n;i++) {if(coprime(i,n)==1) ret++;} return ret;} 
main() 

{ int kk; 


for(kk=1;kk<=N;kk++) if(SumS(kk) ==Euler((kk*(kk+1) /2))) 
printf( “rusult is % d\n” ,kk); 
getch (); } 
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Abstract For orthogonal projective matrix R, i.e., R? = R and RT = R, we say that A is 
generalized Hermitian matrix, if RAR = A*. In this paper, we investigate the least residual 
problem ||AX — B|| = min with given X, B, and associated optimal approximation problem 
in the generalized Hermitian matrix set. The general expressions of the solutions are derived 
by matrix decomposition. 
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81. Introduction 


Some symbols and notations: Let C?"*” be the set of all m xn complex matrices with rank 
r, HC"*” be the set of all n x n Hermitian matrices. Denoted by At, A*, rank(A) the Moore- 
Penrose generalized inverse, conjugate transpose, rank of matrix A, respectively. Moreover, I, 
represents identity matrix of order n, and J = (€n, @n—-1,°+ ,€1), €; € C” is the ith column of 
In. stands for the Frobenius norm. Matrix R € C?’*” is said to be projective (orthogonal 
projective) matrix, if R? = R (R? = Rand R* = R). 

Definition 1.1. If A € C"*”, we say that A is centro-symmetric matrix, if JAJ = A. 


The centro-symmetric matrix has important and practical applications in information the- 


ory, linear system theory and numerical analysis (see [1-2]). As the extension of the centro- 
symmetric matrix, we define the following conception. 

Definition 1.2. For given orthogonal projective matrix R € C?*”, we say that A € C”*” 
is generalized Hermitian matrix, if RAR = A*. Denote the set of all generalized Hermitian 
matrices by GHC”*”. 

In this paper, we discuss two problems as follows: 

Problem I.(Procrustes Problem): Given orthogonal projective matrix R € R"*”, and 

X,Bec"*™, find A € GHC”*” such that 
|| AX — B|| = min. 

Problem II.(Optimal Approximation Problem): Given M € C”*”, find A € Sg such 

that 
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|| M— Al = min [IM — All, 
AE€Sz 

where Sz is the solution set of Problem I. 

Obviously, when M = 0, Problem II is changed into finding the least Frobenius norm 
solution of Problem I. 

Many important results have been achieved about the above problems with different ma- 
trix sets, such as centro-symmetric matrix!!, symmetric matrix!4—®!, R-symmetric matrix!®—7] 
and (R,S)-symmetric matrix!®) set. In this paper, we investigate the above problems in the 


generalized Hermitian matrix set by matrix decomposition. 


§2. Preliminary knowledge 


In this section, we discuss the properties and structures of (orthogonal) projective matrices 
ReEC?P*" and A € GHC™*". 
Denote s = rank(I— R), we know that r+s =n since R? = R. Suppose that pi, po,..., Pr 
and q1,q2,---;@s are the normal orthogonal basis for range R(R) and null space N(R) of R, 
respectively. Let P = (pi, p2,..-,Pr) € C?*" and Q = (q1,42,---59s) € C?”**, then 
P*P=I,, Q*Q=I,, (1) 
RP = P, RQ =0. (2) 
Lemma 2.1.(see [9]) Let matrix A € C”*™ and its full-rank factorization A = F'G, where 
Fecr*",Gecrx™ then A is projective matrix if and only if GF = I,. 


Lemma 2.2. R €C?*” is projective matrix, then 


I, 0 P 
R=(P Q) ap (3) 
0 O Q 
-1 P 
where matrix (p Q) is invertible, and ( P Q) =/ 2 
Q 
If R is orthogonal projective matrix, we have 
I, 0 P* 
R=(P Q) (4) 
0 O Q* 


where |P Q) is unitary matrix. 


Proof. Assume that the full-rank factorization of R is R = PP, we obtain from Lemma 
2.1 and (1) that 
P=P*R, PP =I. (5) 


Similarly, if the full-rank factorization of [— Ris [— R= QQ, we generate 
Q=Q"(- 8), QQ =I, (6) 
since (I — R)? = I — R. Connecting with (1)(2)(5) and (6), we know that (3) holds. The 
equality (4) is obvious since R* = R. 
Lemma 2.3. Given matrices R as in (4) and Ae GHC”*", then 
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A=(P Q) oo _ Vv GeHo"*”, (7) 


Proof. According to Lemma 2.2 and Definition 2.1, it is clear that (7) holds. 
Lemma 2.3 indicates that arbitrary matrix M@ € C”"*” can be written as 


M, Mo P* 
M3 Ms, Q* 


§3. The solutions of Problem I and II 


Given matrices X,B € C"*™, partition 
P* G P* B, 
X= and B= (8) 
Q* X2 Q* By 
where X;, By € C"™*%™ and X29, By € C®*™, 
We need the following two lemmas derived from References [7] and [8], respectively. 
Lemma 3.1. Suppose that matrices X1,B, in (8), then matrix equation A;X1 = B, is 
consistent for A; € HCO", if and only if ByXtX, = B, and X{B, = BY Xi, the general 
solution is 
Ay = A, +U, — OX) KE — Xi XP), 
where A, = (I, — 21) B, x} + (BLX})'(, — 222), VK, € HO™. 
Lemma 3.2. Given matrices X,, B, in (8), then 


min || GX; — B, ||=|| Bid, — Xf X)) | 
GeCrxr 


if and only if G = B, Xj + Ko(I, — X: Xj), V Ke € C"*". 


According to Lemmas 3.1 and 3.2, we obtain 


Lemma 3.3. For the above given matrices X1, Bj, 


min — || A.X1 — By ||=|| Bid, — XPX1) | 
Ai€HC™xr 
if and only if 


XIB AS HAP RBA, (9) 
and the expression of A; is the same as that in Lemma 3.1. 
Proof. | A, Xy = By \|? =|| By = ByX{X + BiX?X = A, Xy \|? 


=|| Bi(i- — XX) |? + || B.xf xX -— AX |l? 


Hence, the least residual can be attained only if By X X, = A,X, which is consistent for 
A; € HC™*" under condition (9) by Lemma 3.3. The proof is completed. 
Based on the previous analysis, Problem I can be solved in the following Theorem. 
Theorem 3.1. Given matrix R as in (4), X,B € C"*™ and the partition (8), then 
min ||AX — Bl? = |B, — Xf X1) ||? + ||Ball?, (10) 
AEGHCrX” 
if and only if (9) holds, at this time 
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0) GH aX RG a) 0 fe 
0 Ui} \Q* 


+ 
a) VK € HOT, 


A (p (11) 


where G = (I, — Se aie +(Bi Xf) (LE - 


Proof. According to the unitary invariance of Frobenius norm, formulas (4) and (7), we 


obtain 
| AX — Bl]? 

2 
= (p Q) le © ae: 

0 of a 

2 

_ G 0 Xy By 

0 0 Xo Bg 

= ||GX1 — By|| ? + || Boll ?. 


Therefore, the problem (10) is equivalent to the following least residual problem 
| GX, — B, || . 


min 
From Lemma 3.3, we know that the ace can be attained if and only if (9), and 
G=C+ (= Xx KU, — XX), 

where kK € HC**” is arbitrary. Submitting G into (7), then (11) holds. 

The following lemma stated from [6]. 

Lemma 3.4. Let LeCcr™, Ac C™!.7T €C™™, and A? = A= A*, Tl? =f, 
then || L — ALT ||= gut, || L — ANT || if and only if A(Z — N)T = 0. 

Let Sz be the solution set of Problem I. We can easily verify from its definition that Sz is a 
closed convex subsets in matrix space C"*” under Frobenius norm. The optimal approximation 
theorem!!®l reveals that Problem II has unique solution, which can be expressed in the next 


theorem. 


Theorem 3.2. Suppose that the given matrix in Problem II is 


M, Mz\ [P* 
M= (P Q) eu 
M3 M4, Q* 

then 

min || M—A|| (12) 

AESzr 
if and only if 

CLR) 4 ae 0 
A=(P Q) ( ina ok 1X7’) . (13) 


0 0 Q* 
where G is the same as that in Theorem 3.1. 
Proof. By using the unitary invariance of Frobenius norm and Theorem 3.1, we obtain 
‘ M, Mp Git =x KG. = ea Oo 
| Mak - 
M3; M4 0 0 
= || (Wh - G) — (h — XA XP) Kp — Xi XP) |? 
+ || Mo ||? + |] Ms ||? + |] Ma ||?, 
then the problem (12) equals to solve the minimum problem 
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We 20) = BR = | 
ania... I hn — G) — ( 1X7) K( 1X7) | 


Moreover, since ||M,||? = || 


a |? + | a \|?, hence the above minimum problem can be 


transformed equivalently as 


M,+M* — 
min || (— PG) Se SOR = Ba) | 
KEHCTX*r 2 
We further deduce from Lemma 3.4 that 
(b- SR = ea =O, eae a ee), (14) 


submitting (14) into (11), we obtain (13). 
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